An extremal model for amorphous media plasticity 
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An extremal model for the plasticity of amorphous materials is studied in a simple two- 
dimensional anti-plane geometry. The steady-state is analyzed through numerical simulations. 
Long-range spatial and temporal correlation in local slip events are shown to develop leading to 
non-trivial and highly anisotropic scaling laws. In particular, the plastic strain is shown to statis- 
tically concentrate over a region which tends to align perpendicular to the displacement gradient. 
By construction, the model can be seen as giving rise to a depinning transition, the threshold of 
which (i.e. the macroscopic yield stress) also reveal scaling properties reflecting the localization of 
the activity. 
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In contrast with crystalline solids, amorphous mate- 
rials display a plasticity which cannot be resumed to 
the motion of well identified defects such as dislocations. 
Consequently, the microscopic description of amorphous 
plasticity is still lacking a consistent framework. Recent 
studies have focussed on the fact that global plas- 
tic deformation is mostly due to local rearrangements. 
Starting from a molecular dynamics study of a bidimen- 
sional Lennard- Jones glass and measurements of the me- 
chanical response under shear stress, Falk and Langer 
pj introduced the notion "Shear Transformation Zones" 
(STZ) having a bistable character to build a mean field 
theory of plastic deformation in an amorphous material. 
Initially drawn by Bulatov and Argon for amor- 

phous solids materials this approach can be extended to 
granular materials or dense suspensions m. In the fol- 
lowing we study a minimal model of plastic deformations 
in disordered media. This model has been proposed in 
the early eighties to describe the fault self-organization 
in seismic regions where mostly the "avalanche" 

properties of this model were studied in order to com- 
pare with the observed power-law distribution of seis- 
mic events (Gutenberg- Richter law). In 0], a quenched 
random distribution of local threshold stress was intro- 
duced which allowed for a mapping onto a random poly- 
mer problem in the limit of vanishing stress drop. The 
analysis we propose focuses on the scaling feature of the 
spatio-temporal organization of the local slip events and 
on the stress-strain characterization. 

We consider a bidimensional material submitted to 
anti plane shear stress. The elastic component of the 
displacement u z {x,y) is thus solution of a Laplace equa- 
tion V 2 u z = 0. The material is discretized on a regular 
lattice, the axis of which are oriented at 45 degrees from 
the displacement gradient direction. The elastic modu- 
lus is assumed to be uniform but we impose a spatially 
frozen disorder for the onset of slip at a local level. After 
a local slip, we renew the threshold stress where the slip 
occured from a random distribution. Bi-periodic bound- 



ary conditions are implemented for the stress and the 
strain, whereas a discontinuity is imposed on the dis- 
placement along the y axis. We choose an extremal dy- 
namics: the external load S is adjusted at each step so 
that only one bond is at the plastication threshold. In 
the spirit of Ref. [0 , this corresponds to a structural re- 
arrangement of a "shear transformation zone" . The lat- 
ter induces a displacement discontinuity along the bond 
and a random modification of the plastication threshold. 
The local stress is redistributed over the material accord- 
ing to the elastic response function. The local stress on 
a bond i is <Tj = S + cr° where £ is the macroscopic 
stress. After a slip Auj at bond j, a° is adjusted to 
a° = a° + AujG(xj_ — Xj). Apart from the periodicity im- 
posed by the boundary conditions, this function G(x) is 
long ranged, decreasing as G(x) oc |x| -2 . Moreover, due 
to the shear boundary condition, the stress redistribution 
is anisotropic. In the longitudinal x direction, bonds are 
loaded while in the transverse y direction they are un- 
loaded. This anisotropy is one distinct feature of the 
model leading to the localization effect to be described 
below. 

The local plastication thresholds j(x,y) are randomly 
chosen according to a uniform distribution between and 
1. Since the elastic modulus is uniform, the stress redis- 
tribution is computed once and for all via a conjugated 
gradient algorithm, for a local slip of unit magnitude. 
This effective Green function is then simply translated to 
the location of the slip event and its amplitude is scaled 
by the slip magnitude (also assumed to be random). 

The behavior of this model is of the pinning/depinning 
type. The plastication criterion of an individual bond can 
be written Yf xt > 7(2, y) — a el (x, y) where a ext is the ex- 
ternal shear stress, 7(2;, y) the local plastication threshold 
and (i el (x,y) the local stress component due to elastic 
stress redistribution from previous plastication events. 
Following an extremal dynamics, we select at time t 
the current "weakest site" (x*,y*) such that a c (t) = 
l{x*,y*)-a el (x*,y*) = mm {x>y )[y(x,y)-a el (x,y)]. The 



1 



maximum over time a* = max t a c (t) corresponds to the 
macroscopic yield stress. From this signal we can recon- 
struct the evolution of the system subjected to a constant 
load: when submitted to an external shear stress lower 
than a* the plate deforms plastically before blocking in 
a jammed state. For values above the yield stress, the 
system flows indefinitely. 

Such pinning systems have been extensively studied 
over the recent years. They have been used to describe 
front motion in a disordered environment in the context 
of wetting 0, magnetic domain walls JlO[ , fluid invasion 
in porous media fTl] ], crack propagation jl2 14 1; or the 
behavior of 'periodic systems', e.g. vortex lattices [ fL5| or 
charge density waves (CDW) [ ^6[ . 

Beyond the yielding transition, this simple model ex- 
hibits another characteristic feature of plasticity: hard- 
ening (i.e. increase of the yield stress with the plas- 
tic strain). In crystalline solids, the hardening behavior 
is due to the entanglement of dislocation loops. In the 
present case, after a first loading, we observe an increase 
of the elastic limit. However the mechanism for this hard- 
ening effect is here of a pure statistical nature. During 
the loading process, the weakest sites are progressively 
decimated. Then the plastic threshold 7(2;, y) is renewed. 
The new threshold is in average larger than the previous 
one. This introduces a systematical bias. When submit- 
ted for the first time to a loading process, the distribution 
of these local plastication thresholds evolves to eventually 
reach a stationary state. On Fig. |], we show the evolu- 
tion of the mean plastication threshold (7(0;, y)) during 
loading. The asymptotic steady distribution seems not 
reached yet on the figure. This hardening effect thus cor- 
responds to a progressive reinforcement of the weakest 
regions. 

On Fig. H (above) we show a map of the cumula- 
tive plastic strain for a system of size 128 x 64 after 
810 5 time steps. We see clearly that the plastic strain 
is non uniform: it is localized within regions elongated 
along the x direction. Focusing on the plastic deforma- 
tion taking place within a finite time window, we show on 
the same figure (below) the appearance of an individual 
localized structure. To characterize quantitatively this 
spatial distribution, we studied the pair correlation func- 
tion of the plastic strain e p (x,y) through Fourier trans- 
forms of the strain map averaged over time. We found 
that the projection of the plastic strain along the x or 
y axis, e\\{x) = (e p (x,y)) y and e±{y) = (e p {x 7 y)) x are 
sclf-affinc profiles with roughness exponents Oi w —0.09 
and £j_ w 0.50. Figure |3| shows the power spectra of e p 
for k x = and k y = 0, where the power-law behaviors 
give directly the cited roughness exponents. 



Such a scaling behavior which characterizes the steady 
state fluctuations of the cumulative strain allows to ana- 
lyze the time evolution of the plastic flow. Let us consider 
two local slip events separated by a time lapse r, and 
record their distance along the x and y direction, noted 
respectively d|| and d±. Averaging over time (at fixed r), 
the probability distribution function of these distances 
p(d, t) reveal two characteristic "correlation lengths" , £11 
and £xj below which p is constant, and above which p 
decays as a power-law with an exponent a\\ or a±_ re- 
spectively. Varying the time lapse r, we observe that 



£11 (t) cxt 1 /*" a(r)ocr 1/2j 



(1) 



Exploiting the self-affine nature of the cumulative plas- 
tic strain, and using a result obtained for other extremal 
models of depinning, we can relate the two dynamic ex- 
ponents to the roughness exponents 1 : 



Z\\ = 1 + C|| z ± = 1 + Ci 



(2) 



The numerical values of the z exponents are consistent 
with these identities. 

The difference in scaling in the a; and y direction can be 
accounted for through a power-law relating both direc- 
tions. Indeed, the correlation lengths are related through 
£j_ c< £jj with = Z\\/z± » 0.65. Moreover, looking at 
the mean value of d± for a prescribed value of d|| also 
reveal the same power-law d± oc djj, with (3 « 0.65. 

Let us focus now on the depinning stress distribution. 
On Fig. H we show the distribution of the local plastica- 
tion stresses, 7(2;, y) — <r el (x, y) (at all sites and all times) 
and of the current plastication stresses, a c (t). The max- 
imum of the latter over time corresponds to the macro- 
scopic yield stress a* . We clearly see that the yield stress 
separates two distinct regions. Stresses larger than a* 
are approximately distributed according to a normal law. 
Stresses lower than a* are however distributed accord- 
ing to a power law of the argument (<r* — a). As above 
suggested the fraction of sites such that a < a* can be 
thought of as a population of potential active sites and 
hence may be interpreted as potential STZ. Let us em- 
phasize that these STZ are not postulated but emerge 
naturally within the model. 

Prior to a large jump in the location of the slip event, 
the lattice has reached a state of strong pinning. Hence, 
following the analysis presented in Ref . |l7j , if we condi- 
tion the statistical distribution of a c (t) by the distance 
to the location of the next slip event, along the x direc- 
tion for instance, Ax, we observe that the larger Ax, the 
narrower the distribution and the closer its mean to the 



1 When £ < 0, an effective value of C, R = should be read 
in this formula 
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yield stress a*. These distributions are shown in Fig. 
H Motivated by the underlying criticality of the depin- 
ning transition, we may anticipate a scaling form of the 
distribution as 



p(a c \Ax) = Ax^ip\ 



)Ax^ 



(3) 



This particular form implies that the standard devi- 
ation of the distribution, 5a c vanishes as (Ace) - "!!, and 
that the mean value of {a* — cr c )(Aa;) is simply propor- 
tional to 6a c (Ax). The first property allows to determine 
and the second gives a simple way to estimate pre- 
cisely a* through a simple linear regression. The same 
procedure applied to Aa; gives a similar result. Using the 
linear dependence of Sa c on a c conditioned to the size 
of the activity jump in both the x and y direction, we 
find numerically a* — 0.517 for a uniform distribution of 
threshold 7 in [0, 1] and a random slip amplitude from 
the same distribution. 

The scaling of the standard deviation of the distribu- 
tion versus the jump size gives a determination of the 
exponents fss 0.68 and v±_ « 0.98. We note that again 
the ratio of these exponents gives the anisotropy scaling 
(3 = v»/v± = 0.69 in good agreement with the previous 
determinations ((3 « 0.65). 

The knowledge of the distribution p(d) cx d _Q n of the 
x-distances between successive active sites allows to ex- 
press the depinning stress distribution close to threshold: 

Q{a* -a c ) = J x v W- a Wil) [{a* - a c )x v W ] dx 

« (a* - <r c r (4) 

where 
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(5) 



The same argument obviously also holds for the y di- 
rection. This latter scaling is also consistent with the 
anisotropy scaling (3 = (an — l)/(a± — 1) « 0.64. 

Despite its extreme simplicity, the model that we pre- 
sented accounts for several features of plasticity in amor- 
phous materials. We could identify a macroscopic yield 
stress. Below this threshold, the material deforms plas- 
tically before blocking in jammed state. Above, it can 
flow indefinitely. This behavior is typical of a pin- 
ning/depinning situation. In the same spirit as the study 
presented in Ref. [lffl, the model exhibits a critical be- 
havior of the plastic stress close to the macroscopic yield 
stress. When submitted for the first time to a shear stress 
we observe a hardening effect. In contrast with crys- 
talline materials, here this effect is of a pure statistical 
nature and corresponds to a progressive reinforcement of 
the weakest regions. In addition to this global hardening 
plastic behavior, the model exhibits a statistical localiza- 
tion. The latter appears via elongated structures in the 



shear direction. However, instead of concentrating onto 
a unique structure (such as in Ref. || ) , the plastic strain 
develops a complex spatio-temporal organization. A sta- 
tistical analysis of these patterns reveals scaling proper- 
ties; scaling exponents are summarized in table [j]. 

Beyond this simplified model, the introduction of ther- 
mal activation in the selection of the site to plastify 
should allow to account for visco-plastic effects. Another 
improvement of such models would consist in including 
both deviatoric and volumetric strain, the latter coupling 
being characteristic of irreversible deformation in amor- 
phous solids. 
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C|| = -0.09 ±0.05 


Ci = 0.50 ±0.05 


a|| = 1.61 ±0.05 


a± = 1.96 ±0.05 


v\\ = 0.68 ±0.05 


v± = 0.98 ± 0.05 


/i = 0.00 ±0.05 


/3 = 0.65 ±0.05 
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TABLE I. Table of scaling exponents 



4 



10' 



10" 



10 



15 



20 



<£ p> 



FIG. 1. Evolution of the mean plastication threshold 
(7(3;, y)) during the transient regime of a first loading process. 
The increase of (7(1,2/)} can be interpreted as a hardening ef- 
fect. 
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FIG. 3. Power spectra of the plastic strain e p for k x = 
(circles) and k y — (squares). The lines indicate power law 
behaviors corresponding to roughness exponents <ji = —0.09 
and Cl = 0.50 




FIG. 2. Map of the relative displacement field obtained on 
a system 128 x 64 after 800 000 times steps (above) . The dif- 
fuse localization corresponds to the succesive development of 
anisotropic structures elongated in the longitidinal direction. 
Focusing on a finite time window (450 time steps, below) al- 
lows to reveal an individual structure. 



FIG. 4. Distribution of the individual site plastication 
thresholds Ui (line) (all sites, all times) and of the current 
plate plastication thresholds a c (dotted line) (active site, all 
times). The vertical line indicates the position of the yield 
stress a* . 
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FIG. 5. Distribution of depinning stress (bold) and contri- 




butions conditionned by the distances (2,4,8,16,32) between 
consecutive active sites. The tail of the distribution corre- 
sponds to very short jumps and is very sensitive to the de- 
tails of the random threshold distribution. The contributions 
obtained for increasing distances between consecutive active 
sites present the same trend: the larger the jump, the closer 
the mean force to the threshold and the narrower the distribu- 
tion. After rescaling (inset) these distributions collapse onto 
a single master curve. 
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